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^ ■ ABSTRACT 

^-H I The relevance of orbital eccentricity in the detection of gravitational radiation from 

(steady-state) binary stars is emphasized. Computationally effective (fast and accu- 

I rate) tools for constructing gravitational wave templates from binary stars with any 

■ orbital eccentricity are introduced, including tight estimation criteria of the pertinent 

' truncation and approximation errors. 

o : 

in ■ 

o : 

(3 1 INTRODUCTION 



o 



Gravitational wave detection experiments in space, including satellite Doppler- Tracking (Bertotti and less, 1999) and LISA 



(tittp:/ /lisa.jpl.nasa.gov), will hopefully open a window on the low- frequency part of the gravitational wave (henceforth GW) 



I spectrum of cosmic origin. In these frequency bands, binary stars are among the most promising continuous detectable source. 



A substantial fraction of binaries are expected to have orbits with non negligible eccentricity (Barone et al., 1988; Hils et 
al., 1992; Pierro and Pinto, 1996c) resulting into the emission of several harmonics of the fundamental orbital frequency. The 
^ importance of this fact from the standpoint of signal detection and estimation has been already noted. 

, For coalescing binaries, Pierro and Pinto (1996b) and Martel and Poisson (1999) pointed out that neglecting residual 

5_j ' (albeit very small) orbital eccentricities may seriously deteriorate matched-filter detection performance. Their results, obtained 
in the frame of the simplest (newtonian) Peters Mathews (henceforth PM) model (Peters and Mathews, 1963; Peters 1964; 
Pierro and Pinto 1996c), support the qualitative conclusion that residual orbital eccentricities cannot be bona fide disregarded 
in building templates for matched-filter detection of gravitational wave chirps from inspiraling binaries^ 

For steady-state binaries with non-zero orbital eccentricity, on the other hand, using circular-orbit waveform templates, 
i.e. neglecting higher order harmonics, implies a potentially large loss of signal-to-noise ratio (henceforth SNR), leading to 
significantly worse detector's performance, as will be shown in the sequel. 
The main goals of the present paper are: 

i) to provide some quantitative hint for validating the applicability of the simple PM model to steady-state binaries; 

ii) to gauge the loss in SNR due to the simple circular-orbit assumption and, more generally, to set some criteria for spectral 
waveform truncation; 

iii) to introduce efficient (accurate and fast) computational tools for constructing gravitational waveform templates for 
(steady-state) binary sources with any orbital eccentricity. 

The paper is organized as follows. In Sect. 2 we introduce some (dimensionless) parameters whereby the applicability 
of the PM model to specific sources can be assessed. In Sect.s 3a and 3b we review the GW spectra and waveforms in the 
frame of the PM model. In Sect.s and ^6 we show how to evaluate the total harmonic distortion due to spectral waveform 



* The effect of a residual (tiny) orbital eccentricity on the radiation emitted from an inspiraling binary system was also considered in 
(Moreno et al., 1994), with special emphasis on the possible relevance of periastron advance. The results in (Moreno et al., 1994) are 
unfortunately affected by several errors and misprints. 
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truncation, and introduce a modified Carlini Meissel expansion tool for fast and accurate GW harmonics computation. The 
results in this section can be readily extended, in principle, to higher-order post-newtonian (henceforth PN) models. As an 
application, in Sect. 5 we apply our formalism to some paradigm eccentrical binary sources. Conclusions follow under Sect. 6. 
Technical developments are collected in Appendix A to C. 



2 STEADY-STATE BINARIES: THE PETERS-MATHEWS MODEL 



The PM model for gravitational wave emission from binary systems in a Keplerian orbit was introduced in the sixties (Peters 
and Mathews, 1963, Peters 1964), and recently re-examined (Pierro and Pinto, 1996a) . It relies on the following main 
assumptions: i) point mass, ii) weak field, iii) slow motion, and iv) adiabatic evolution (negligible change of the orbital 
parameters over each orbit). These conditions can be checked in terms of the following inequalities (Pierro and Pinto, 1996a): 
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where as already stated x = cT/ivrg, T being the orbital period, rg = 2G{Mi + M2)/(? is the source gravitational radius. 
Ml, 2 are the companion masses, A = \M\ — M2\/(M\ + M2) and e is the eccentricity. 

Tidal effects could be neglected provided neither companion star fills its Roche lobe. Following Eggleton (1983), this 
translates into: 
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where A is the ratio between the physical and gravitational companion radiusfj. 

For most steady-state binary systems, i.e. long before coalescence, ,^1 to ^3 above are fairly small (see e.g. Sect. 5), and 
the PM model turns out to be perfectly adequate. 



3A STEADY-STATE BINARIES: SPECTRA 

According to the PM model, the GW power C^'^ radiated at the n*'' harmonic of the orbital frequency by a steady binary 
source can be conveniently cast into the following universal form (Barone et al., 1988): 

cl'"' = |gx-^°^^(l- A^)'G_(e)ff+'^(n,e) (5) 
be" 

where the superscripts -f-, x refer to the fundamental GW polarization states. The spectral power distribution is embodied in 
the universal dimensionless functions g~^'^ (n, e) shown in fig.s 1.1-1.10 for e = 0(0.1)0.9. For circular orbits (e = 0) only the 
second harmonic is emitted. The function Gmax{e) plotted in fig. 2 is the ratio between the total luminosity (sum over both 
polarizations) of the brightest GW spectral line, and the total luminosity of a circular-orbit binary having the same x ^-id ^■ 
The brightest spectral line is the Nmax-t\i harmonic of the orbital frequency, where Nmax is a function of e only, displayed in 
fig. 3. 

It is seen that for non circular orbits, several spectral lines with comparable intensities are emitted. Thus, use of the 
circular orbit waveform templates implies a potentially sizeable loss in the available signal power and hence in the SNR, which 
can spoil the detector's performance. 

3b WAVEFORMS 

The far-field metric deviation (TT gauge) in the PM model is^: 

^ Typical values of A range from 10^ for white dwarfs down to 3 for hadronic stars. Further departures from the standard model are 
expected due to the possible occurrence of mass-transfer phenomena, which would be present in closely-orbiting classical stars, as well 
as in binaries where one companion is an accreting collapsed object. 

i The GW field can also be obtained by inverting the Keplerian integral of motion relating time to the true anomaly, and exploiting 
the simple dependance of the radiated waveforms on this latter (Wahlquist, 1987). The referred procedure is purely numerical and, to 
the best of our knowledge, its generalization to higher PN order models is not immediate. 
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COS '(^ 

hx = — ^ [2ha:y cos2(fi — {hxx — hyy ) sin 2(p] , (6) 
V 2 

1 f3 + cos2i9 . „ , , „ , 1 — cos2i?,, , ,1 
h+ = — I [2hxy sm2(p + (h^x - hyy) cos2<p\ {ha:^: + hyyjj , (7) 

where the coordinates i9 and i/p specify the direction of the observer in a spherical polar system where the orbit lies in the 
equatorial plane and the binary center of mass is at the origin. 

The metric components in can be expanded into Fourier series under the adiabatic assumption that the orbital 

parameters do not change appreciably over each orbit. Hence|^ 

^ 2 

h^y = ^ hi"J sin [n^tj y (8) 

n = l 

^ 2 

K±y = h^"2y COS (n^t^ , (9) 
where hx±y is a shorthand for hxx ± hyy (see Appendix A), 

hi"J = hon{l - e)^^'^ [Jn-2(ne) + J„+2(ne) - 2J„(ne)] , (10) 
= 2/iow{J„_2(ne) - J„+2(ne) - 2e [J„_i(ne) - J„+i(ne)] + (2/n) J„(ne)} , (11) 
h^^ly = -4/ioJn(ne), (12) 
dPI 



au' 
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/lo = 1 T-pT- (13) 

47rr X 

For circular orbits one has simply: 

'^i-^. = 2e;=4;jo5„2, fti;'.=0, (14) 

where is the Kronecker symbol. 

For steady state binaries the (Robertson) periastron advance|||] does not produce sensible effects on the waveforms, and is 
thus deliberately ignored. Inclusion of the periastron advance amounts to splitting each GW spectral line into a doublet at 
~ (27r/r)(l ± (3x~^^^), which cannot be resolved unless the signal is Fourier-transformed over a timespan ~ -yP-I^T sec. This 
time is, e.g., ~ 5 ■ 10^ years and ~ 2.8 • lO'^ years for P5'i?1534 + 12 and PSRVdVi + 16, respectively. 



4A SPECTRAL TRUNCATION AND APPROXIMATION ERROR 

In order to discuss the effect of spectral truncation of (^) and ^ on the available SNR it is convenient to introduce the total 
harmonic distortion (henceforth THD): 
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\h- h 
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i,{")l2 



(15) 



where h, h represent the exact and approximate values of the metric tensor, /i'"', h are the Fourier coefficients of h, h, 
respectively, and the L^-norms are computed by taking the time average over one orbital period of the square of the argument, 
within the spirit of the adiabatic approximation. If only A'^r harmonics are included, then 



§ The unknown irrelevant phase at t = has been set to zero. 

^ Note that for n = 1 eq.s ( px| ) and (jll|) contain Bessel functions of order —1, for which J—i(x) = —Ji(x). 

II The relativistic periastron advance was heuristically (i.e., inconsistenlty, from the post-newtonian expansion view point) included in 
(Moreno et al., 1995). 
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0, n > Nt 
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(16) 



It is readily recognized that THD^ represents the fraction of signal power which is lost as an effect of truncation^. 

In the most general case, where besides spectral truncation, the Fourier coefficients are computed in approximate form 
(as e.g. in the next subsection), one has: 

/ ATT , , NT , , \ 1/2 

2 2J h -}2 I h I' 

(17) 
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)|2 



The harmonic distortions TY{T)x±y, THD^^ due to the spectral truncation of (|^, (|^ can be computed for any given A'r using 
Kapteyn's theory (Watson, 1966, ch. 17) to evaluate in closed form the infinite sums in (^^. After some lengthy but simple 
algebra, one obtains (see Appendix B): 
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(20) 



The corresponding harmonic distortions for the TT metric components h+, can be conveniently written as follows 



THDx = i^[iTHDly\\hxyf cos'' 2ip + THDl^y\\hx-yf siii' 2ip] ■ [4||/i,j;f cos' 2^9 + sin' 2.^] '| , 



(21) 



and: 

THD+ = I (3+cos2i?)^ {ATHDljh^yfsm''2ip+THDl_Jh^^yfcos''2^)+{l-cos2^fTHDl+y\\h,^^^^^ + 
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where (•, •) is the scalar product in In order to evaluate (|2^) the further infinite sum: 

OO 

= E'^i-^it - -8(l-e')-^/' {1+ (l-2e-') [l - (1 - e')^/'] } , 



(22) 



(23) 



is needed, which is also readily obtained as explained in Appendix B. 

The harmonic distortions ( [ill ) ^^i^ ^^^n be sensible even at very low eccentricities (e < .1). Expanding ( |2l| ) and ( |22| ) 
to lowest order in e yields: 

(24) 
(25) 



^4(1 - cos 2o3)2 + 12 1 - cos 2(;9)(3 + cos 2(^9) cos 2i9 + 90(3 + cos 2(z))2 , 
+ 4(3 + cos2(^) ^ ^ 



The THD is closely related to the fitting factor FF (Apostolatos, 1996) between the exact and spectral-truncated (template) waveform. 
From the very definitions one gets: 

THD 2 



FF ^ I 



0( THD ). 
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The above simple expressions are fairly accurate for e < .1, as seen, e.g., from fig. 4, where the angular averages of the 
approximate and exact harmonic distortion are drawn, and seen to be almost indistinguishable and non-negligible. The 
(i?, (ys)— dependent factor in ( |25| ) is plotted in fig. 5. Its average value over the sphere is exactly equal to the ^)— independent 
factor in (p^. 

The obvious question is how many terms should be included in and so as to keep both THD+ and THDx below 
some specified level, for any 

To answer this question one may resort to the following inequalities: 

maxTHDx < max{THD^y,THD^-y), ma.xTHD+ < max(THD^y,THD:,^y,THD^+y) max[Q(i9, 93, e)]. (26) 

(i9,<^) (i5,c/)) 

where: 

Q(i9,(^,e) = | (3+cos2i9)^ (4||/!.,y f sin^ 2v3+||/i:.-y |1^ cos^ 2v5) +(l-cos 2i9)^ \\h^+yf + 



+ 2|(l-cos2i9)(3+cos2i9)cos2v3| 



3 +COS 2i9) ^ (4|| Kyf sin^ 2ip+\\h:,-y f cos^ 2(p) + 



(27) 



+ (l-cos2i9)^r7i"£)^+y|l/i^+j; 11^ + 2 (l-cos2i9)(3 + cos2i9) cos 2v3 (h^-y,h^+y^ 

The first of ^ follows immediately from (|l|); the second one is obtained from ( [22|) using Schwartz inequality. The supremum 
of the function Q(i),ip, e) occurs at i9 = 7r/2, if — rmv, V e, where Q{tt/2, rmv, e) ~ 1.5 (see fig. 6). 

The truncation orders required to keep THDx,+ < 0.01, deduced from d26h are collected in Table-I. 



eo Nt 



.1 


4 


.2 


6 


.3 


8 


.4 


11 


.5 


15 


.6 


22 


.7 


36 


.8 
.9 


68 
206 



Table I - Truncation orders needed to keep THD+^x < -l- 



4A A GENERALIZED CARLINI-MEISSEL FORMULA 



A key issue for an efficient computation of waveform-templates based on (|l^), ( pT[ ) and ( p^ ) involves clever evaluation of terms 
like: 



J„(ne), Jn±i(ne), Jn±2(ne). 



(28) 



It is well known that, in general, whenever the argument and the order are close (here, in fact they are proportional through the 
orbital eccentricity e), numerical computation of Bessel functions either by series summation (Abramowitz and Stegun, 1968, 
ch. X), or by (re-normalized, downward) recurrence (Press et al, 1992, Sect. 6.5) is inefficient. As a convenient alternative, 
we suggest the following generalization of the well-known (see Watson, 1976, ch. XVII) Carlini-Meissel (henceforth CM) 
expansion: 



J„±fc(ne) ~ 4'^")(ne)*±fc(n,e), 

where (see Appendix C for the detailed deduction) 
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Using ( pi] ) to evaluate the Fourier coefEcients h does not significantly spoil the accuracy of the waveforms. Indeed, 
spectral truncation according to Table-I still yields THD values below 0.1. 



5 PROTOTYPE SOURCES 

As an application of the above, we wrote a code for waveform template construction, and used it to compute the waveforms 
for several prototype sources. Taylor et al. (1993) provide data for 24 binary pulsars. In Table II we quote PSR 1913+16 
and PSR 1534+12, as possible paradigm sources for space detectors, being respectively the most popular and closest known 
binary pulsars. 



Binary 


1534+12 


1913+16 


Right ascension 51950 


15:34:47.686 


19:13:12.46769 


Declination B1950 


+ 12:05:45.23 


+ 16:01:08.0323 


Orbital inclination i [degrees] 


74 


45 


Distance [kpc] 


0.68 


7.13 


Projected semimajor axis ai sin i[light ■ s] 


3.729468 


2.3417592 


Eccentricity e 


0.2736779 


0.6171308 


Orbital period Pb [d] 


0.4207372998 


0.322997462736 


Companion masses [Mq] 


1.34, 1.34 


1.42, 1.41 


(10-^) 


3.4849 


4.3102 


6 (10-3) 


1.3200 


1.4680 


Ca (10-'") 


7.6549 


13.023 


A 





3.5336 X 10-3 


X 


434777882.4767 


316085232.7313 


ho (10-23) 


16.518 


2.0575 



Table II - Paradigm compact binary sources 

The gravitational waveforms at -3 = ip = 0,45,90 deg, computed using 8 harmonics for PSR1534+12 and 22 harmonics for 
PSR1913+16 (consistent with Table-I) are displayed in fig.s 7.1-7.15, and fig.s 8.1-8.15, respectively. By comparison, the 
waveforms corresponding to e = are also drawn. 



6 CONCLUSIONS 

The main results in this paper can be summarized as follows. Orbital eccentricity should not be neglected in detecting gravi- 
tational waves from steady-state binaries, for which the simple Peters Mathews model has been shown to be accurate enough. 
GW spectral truncation criteria have been discussed, and computationally efficient tools/techniques have been introduced 
for constructing rehable templates. We stress that the above tools/techniques could be readily extended, to higher order FN 
models with relative ease. 
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APPENDIX A: RELEVANT TO EQ.S (10) TO (16). 



In the weak-field slow-motion approximation, the cartesian far-field harmonic-gauge metric tensor deviation components in 
(8), (9) are simply related to the source quadrupole tensor lij through: 



_ 2G d^hv , 2G d^hx , 2G d^l 

where 



_ - , _ U, ±XX , _ U. ±yy 



Ixx = HP^ cos'^{4>), lyy = np'^ 8m^{<j)), /a,y = /up^cos(^)sin(0), {A2) 
p being the companion star separation, e the eccentricity, the true anomaly and n the reduced mass. 

The relevant terms of the (reduced) quadrupole moment can be conveneintly rewritten: 

Ixx = tia £, , lyy = fia rj , Ixy = p-a ^r], {A3) 

where a is the orbit semimajor axis, 

^'{^)- "-{^)- 

Then, using the well known Keplerian equations (see, e.g., Watson, 1976, ch. XVII) 

p cos Cj) „ p Sin0 2n1/2 . / 

^=cosS-e, - = {l-e)' smE, (A5) 

where E is the eccentric anomaly, and the relation between the latter and the mean anomaly M, 

M = = E-esinE, {A6) 
one can expand ,rf and ^r; into Fourier series of argument M, taking properly into account their parities, viz.: 

oo 

= Y + $^7nCos(nM), (^7) 

n=l 



2 

71 = 1 
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iV = ^Vn sin (nM) . (^9) 



The relevant Fourier coefficients are readily found. Hence, using (^45) and {A6): 



2 r 

7„ = — / sm[n{E - esmE)]{ain2E - 2esmE) dE, (AW) 
"^^ Jo 



5n = — —{l-e^)[ sm[n(E -esmE)]sm2E dE, (All) 
Jo 

Tin = — (1 - e^)^/^ / cos[n(S - esmE)](cos2E - ecosE) dE. (A12) 
Jo 

Upon repeated use of trivial trigonometric identities, and in view of the integral definition of the Bessel function of the 1st 
kind, 



L{a) = 

^ Jo 



cos [i^a; — a sin a;] da, (^13) 



the Fourier coefficients (^10) to (^412) can be written: 



1 2e 
7n = - [Jn-2(ne) - J„+2(ne)] - — [Jn-i(ne) - J„+i(ne)] , (^14) 

1 2 

6„ = --{l-e ) [Jn-2(ne) - J„+2{ne)] , (^15) 

r?„ = -(1 - e^Y^'^ \Jn-2ine) + J„+2(ne) - 2J„(ne)l . (^16) 
n 

Using (^14) to (A1&) and [Al) to (AQ) in (Al) to (A3) gives equations (10) to (16). 



APPENDIX B: RELEVANT TO EQUATIONS (22) - (26) 

In order to establish eq.s (22) to (26) one may repeatedly use the recurrency formula: 

Jn±l{z)='^JJz)±j'„(z), {Bl) 

so as to reduce the sought series to combinations of the following (generalized) Kapteyn's expansions of the second kind: 

oo 

Y,n[j'n{ne)]\ (52) 

■n=l 

oo 

Y,[jUne)\\ (53) 

n=l 

n=l 

oo 

^nV^(ne), (B5) 

oo 

^J^(ne), (B6) 

n=l 

oo 

^^nJu{ne)Jn{ne). {B7) 

n=l 

These latter can be summed as follows. Prom the Fourier analysis of Kepler motion, the following equations are readily 
established (see, e.g., Watson, 1966), ch. 17.2 : 
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cosS = -I + 2 V ^^^2^ cos nM, (B8) 
2 ^ — ' n 



n = l 



sm S = - y sin nM, (59) 

n = l 



dM 

Differentiating eq. (-B8) w.r.t. M , and using (BIO), one gets 

siniS 



^ = (l-ecos£)-\ (BIO) 



1 — e cos E 
cos E — e 



2^X(ne)sinnM, (Bll) 

n = l 
oc 

2 ^njl(ne) cos nM. (B12) 



^ J„(ne)cos(nM), (B13) 



(l-ecos£;)3 

Similarly, from (-B9): 

cosg _ 2 
1 — e cos E e 

n = l 
oo 

— — TXT = -V nJ„(ne)sin(nM), (B14) 

(1 — ecosi5)3 e 

n — 1 

where E is the eccentric anomaly, M the mean anomaly, and e the eccentricity. The following procedure can be then applied 
to eq.s (-B8) and (i311)-(i314): i) squaring; ti) taking the average in Ai over (0,27r), using again eq. (BIO); in) using the well 
known (Euler) transformations: 

cosE = {z + z'^)/2, sin E = -i{z - z'^)/2, dE = -iz'^dz, 

so as to express the sought series as contour integrals on |zj = 1 of rational functions of z, which are trivially computed in 
terms of residues. 

As an example, applying the above procedure to eq. (B13), one gets: 



EJ^ (ne) — — — cos E _ / — {1 + z ) ^ _ 

_^ 27r 1-ecosB 2Tvi J^^^^^ z^l^z ~ 2e{l + z^)] 



y Res 

|zil<l 



2\2 



(1 + ^^) 



22[42- 26(1+^2 



(B15) 



The integrand function on the r.h.s. of (B15) has a double pole at 2 = and two simple ones at 2 = (2e)~^[l =F (1 - e^)^/^]. 
Only two poles above fall within \z\ < 1, and (B15) gives: 

00 

n=l 

in agreement with Watson, ch. 17.6, eq. (2). Similarly, starting from (B8), (Bll), (B13) and (B14) one gets, respectively 

00 



tt Note that equation (3) in Watson ch. 17.6, is in error, as seen by comparison with (520), and by direct numerical check. For this 
(erroneous) result Watson quotes (Schott, 1912). The same error appears in (Prudnikov et al., 1986, sect. 5.7.31). 
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The series (-B7) can be summed by differentiating both sides of eq. (-B16) w.r.t. e. Hence: 



Y,nMne)Une)^ \ (B21) 

— ■ 4 ( 1 — I ' 



APPENDIX C: GENERALIZED CARLINI-MEISSEL EXPANSIONS 



To obtain the generalized Carlini Meissel expansion for J„±k{ne) we start from Bessel equation for J„±k{ne): 

d?J„±k{ne) ldJn±k{ne) 



+ 



e de 



+ 



^2 {n±kY 



and let0: 



Jnine) = Jn±kine) = 



(n±fc)! 



exp 



Jn±k{ne) = 0, 



Un±k{z) 



On letting eq.s (C2) into (CI), we get: 



2 -1 2 i 'i^) 

Unifc + u„^k + e iinifc + n 5 = , 

then, following Carlini and Meissel, we assume that the following asymptotic representation for Un±k, {k — 0, 1, 2) holds 

K±kiz) 



Un±kiz) 



n 



Substituting (C4) into (C3), and equating like powers of n (as required by consistency), we get: 



life 



2e2 u+ 



Hence: 



±fc 



2(1 -z2) z./T^ 



(1-^2)5/2-^ (1-^2)2' 2(l-z2)5/2- 

Carrying out the integrations in (C2), and taking into account that Jri±fc(0) = Sn±k,o we get: 



^n±fc(2) dz = y/l-e^ + log 



Wn±fc(2) dz = - ^ log (1 - e^) ± fc 



1 + 



+ c: 



1 + Vl-e2 

tt This formula is suggested by the well-known McLaurin expansions of Bessel functions. 



+ C„±k, 



(CI) 

(C2) 

(C3) 

(C4) 

(C5) 
(C6) 

(C7) 

(C8) 
(C9) 

(CIO) 
(Cll) 
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Plugging the last three eq.s into eq. (C2) we obtain: 



(C12) 



Jn±fc(ne) = ■Pf7TT(l + ^1 - e2)-("±'=) (1 - e')-^/^ • exp i n^l - + n 



(n ± k) 



-ie^ - 2 
24(1 - e2):V2' 



2(1 -e2) 2(l-e2)3/2 

The unknown integration constants can be found by enforcing the following obvious asymptotic equality, valid for all n: 



(C13) 



Jn±k (e — » 0) 



(ne/2) 



(n±/t) 



(n±fc)! 



(C714) 



Hence: 



' 1 + C+, . = 0, 

' n±k ' 

[ -l/12Tfc/2-A;V2 + C-±,=0. 



(C15) 



Hence, from (C13): 



Jn±k(ne) 



(n±fc)! 



-(n±k) 



(1 — e ) ' • exp < n 



V'l - e2 - 1 



+ n" 



-3e^ - 2 
24(l-e2)3/2' 



■ A;^ ttt:^ :rT^ + t:^ ± ^ + ■ 



(C16) 



2(1 -e2) 2(l-e2)3/2 ' 12 2 2 

The r.h.s. of eq. (C16) above will be henceforth denoted as J^-^(ne), and can be more conveniently written as in (29) to (31). 
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Fig. 1.1 - g + x- (n,.0), relevant to eq. (5). 
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Fig. 7.1 - PSRl 534+12 - a waveform gallery. 
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Fig. 1.2 - g + x- (n,.l), relevant to eq. (5). 
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Fig. 2 - Gniax(e), relevant to eq. (5). 
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Fig. 8.2 - PSR1913+16 - a waveform gallery. 
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Fig. 1.3 - g+^x- (n,.2), relevant to eq. (5). 
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Fig. 3 -Nmax (e), relevant to eq. (5). 
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Fig. 1.4 - g + x- (n,.3), relevant to eq. (5). 
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Fig. 8.4 - PSR1913+16 - a waveform gallery. 
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Fig. 1.5 - g +^x- (n,.4), relevant to eq. (5). 
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Fig 5 - The angle-dependent factor F in eq. (25). 
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Fig. 1.6 - g +^x- (n,.5), relevant to eq. (5). 
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Fig 6 - The function Q((?,^,e), eq. (27). 





i3 


= 45(^=90 









12 3 




Fig. 8.6 - PSR1913+16 - a waveform gallery. 
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Fig. 1.7 - g +^x- (n,.6), relevant to eq. (5). 
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Fig. 8.7- PSR1913+16 - a waveform gallery. 
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Fig. 1.8 - g+^x- (n,.7), relevant to eq. (5). 
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Fig. 1.9 - g +^x- (n,.8), relevant to eq. (5). 
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Fig. 8.9 - PSR1913+16 - a waveform gallery. 




Fig. 1.10 - g+^x- (n,.9), relevant to eq. (5). 
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Fig. 7.10 - PSR1534+12 - a waveform gallery. 
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Fig. 8.10 - PSR1913+16 - a waveform gallery. 
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Fig. 7.11 - PSR1534+12 - a waveform gallery. 
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Fig. 8.11 - PSR1913+16 - a waveform gallery. 
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Fig. 8.12 - PSR1913+16 - a waveform gallery. 
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Fig. 8.14 - PSR1913+16 - a waveform gallery. 



